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How did the Sun affect the climate when life evolved on the Earth? - A 
case study on the young solar twin t^} Ceti 
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Using k} Ceti as a proxy for the young Sun we show that not only was the young Sun much more effective in protecting 
the Earth environment from galactic cosmic rays than the present day Sun; it also had flare and corona mass ejection rates 
up to three orders of magnitude larger than the present day Sun. The reduction in the galactic cosmic ray influx caused by 
the young Sun's enhanced shielding capability has been suggested as a solution to what is known as the faint young Sun 
paradox, i.e. the fact that the luminosity of the young Sun was only around 75% of its present value when life started to 
evolve on our planet around four billion years ago. This suggestion relies on the hypothesis that the changing solar activity 
results in a changing influx of galactic cosmic rays to the Earth, which results in a changing low-altitude cloud coverage 
and thus a changing climate. Here we show how the larger corona mass ejection rates of the young Sun would have had 
an effect on the climate with a magnitude similar to the enhanced shielding capability of the young Sun. 
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1 Introduction 

Stellar evolution models predict that when life evolved on 
the Earth more the 4 billion years ago just after the heavy 
bombardment, the conditions in the inner part of the solar 
system were significantly different compared to what we 
know today. E.g. it was noted akeady in 1972 by Sagan 
and Mullen that the luminosity of the young Sun was only 
around 75% of its present value, which would result in 
freezing temperatures on the Earth - assuming a radiation 
budget similar to what we have today (Sagan & Mullen 
1972, see also Gough 1981). This problem has since been 
known as the faint young Sun paradox. On the other hand 
stellar evolution models predict that the activity of the Sun 
related to the chromosphere and corona should be much 
stronger than what we know today. 

A number of different solutions have been proposed 
to the faint young Sun paradox. Sagan & Mullen (1972) 
suggested that elevated levels of CO2 could have main- 
tained surface temperatures above freezing, but the CO2 
level needed to raise the temperature might be so high that 
it would be in conflict with geochemical records (Rye et 
al. 1995). Potentially other greenhouse gasses like ammo- 
nia could also help to raise the temperature, but on the other 
hand it is questionable if large amounts of ammonia could 
be maintained in the Earth atmosphere at a time where the 
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Sun's UV radiation was up to 10 times larger than today 
(Sagan & Chyba 1972). 

Another solution that was originally proposed by 
Graedel (1991) suggests that the Sun has experienced a 
mass loss of 5-10% over its main-sequence life. There- 
fore the young Sun would have been a bit more massive 
than today's Sun and thus brighter than what we would ex- 
pect without considering mass loss. Unfortunately mass loss 
rates as high as 10% contradicts both solar evolution mod- 
els calibrated using helioseismology (Guzik & Cox 1995) 
and measurements of stellar winds around solar-type stars 
(Gaidos et al. 2000; Wood et al. 2002). 

Recently attempts to solve the faint young Sun paradox 
have been based on atmosphere models including not only 
an increased greenhouse effect, but also a reduced albedo 
(von Paris et al. 2008; Kitzmann et al. 2010). These at- 
tempts are in line with the hypothesis put forward by Svens- 
mark & Friis-Christensen (1997): that galactic cosmic rays 
(GCRs) modulate the amount of aerosols and clouds in the 
lower part of the Earth's atmosphere. In other words, in or- 
der to understand the evolution of the Earth's climate now 
and back then, it is important to understand not only the ef- 
fect from high-altitude clouds (through the greenhouse ef- 
fect), but also the effect from low-altitude clouds (through 
the albedo effect). By including both effects in there modes 
it was suggested by Shaviv (2003) and Svensmark (2003, 
2006) that the young active Sun's increased ability to pro- 
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Table 1 Stellar parameters for Ceti (from Gaidos & 
Gonzalez 2002) and the Sun (from Christensen-Dalsgaard 
et al. 1996) 



Name 


Type 


Teff [K] 


logs 


[Fe/H] 


Ceti 


G5 V 


5747 (49) 


4.53 (0.06) 


0.11 (0.04) 


Sun 


G2V 


5778 


4.44 


0.00 



tect the Earth from GCRs could cause higher temperatures 
on the Earth. 

The hypothesis is that as the Sun was much more active 
when Ufe started to evolve on the Earth (which is reflected 
in its higher rotation rate and higher level of UV and X- 
ray emission) it would have been much more efficient in 
shielding us from the GCRs, which would have resulted in 
smaller amounts of aerosols and clouds in the lower part 
of the Earth's atmosphere and thus higher temperatures on 
the Earth. This hypothesis was recently strengthened by ob- 
servations of direct evidence of a relation between GCRs, 
aerosols and clouds on short time scales as observed during 
major Forbush decreases (Svensmark, Bondo & Svensmark 
2009). 

Though it was noted by Svensmark, Bondo & Svens- 
mark (2009) that large Forbush decreases today are too rare 
to have any significant effect on the Earth's climate, this 
might not have been the case 4 bilhon years ago when life 
started to evolve on the Earth as CMEs are expected to have 
been much more common on the early Sun. We therefore 
analyze if the reduction in the influx of GCRs originating 
from Forbush decreases could be significant compared to 
the reduction in the influx of GCRs originating from the 
more effective shielding capacity of Sun at the time life 
evolved on the Earth. 

We undertake this analysis through a case study of the 
young solar twin Ceti . With an age of around 700 mil- 
lion years Ceti mimics the Sun at that time. The age es- 
timate of Ceti is based on the rapid rotation of Ceti 
with a period of 8.6 days (Rucinski et al. 2004) and a re- 
sulting large activity level (Baliunas et al. 1995), but the 
uncertainties of such a simple scahng relation for the ages 
are of course huge. Other fundamental stellar parameters of 

Ceti such as effective temperature, surface gravity and 
metaUicity come so close to solar values that Ceti qual- 
ify as a solar analogue (see Table 1). By using the well stud- 
ied Ceti as a case study instead of (simple) scahng laws 
between stellar activity and age [as done by Shaviv (2003) 
and Svensmark (2003, 2006)] we can base our analysis on 
actual measurements of steUar activity rather than estimates 
based on physical assumptions. 

2 The effect of more effective shielding 

The temperature response to a change in the GCR influx is 
given by (Shaviv, 2003): 



where e„icoti is the GCR flux reaching the troposphere 
around an Earth-like planet around Ceti and Eq is the 
GCR flux reaching the troposphere around the Earth. D and 
q are constants (Z) ^ 10 K and q ^ 0.5). 

The GCR influx can be found by solving the spherically- 
synnmetric transport equation for the stellar modulation 
of cosmic rays reaching the planet's troposphere (Perko, 
1987): 



dU VP dU 

+ 



0, 



(2) 



dr 3k dP 

where U is the cosmic ray distribution function, r is the he- 
liocentric radial distance, P is the particle rigidity, V is the 
solar wind speed and k is the diffusion coefficient for radial 
propagation. The spherically-symmetric transport equation 
can be solved for GCRs with energies larger than a few GeV 
using the force-field approximation. In the force-field ap- 
proximation the kinetic energies of the GCRs at the planet 
E is given as: 

E - EisM - (3) 

where ii^isM is the kinetic energy of the GCRs in the in- 
terstellar medium at the astrospheric boundary and $ is the 
modulation strength (Perko, 1987): 



$ = — . 

3k 



(4) 



If we assume that the stellar wind speed and the diffusion 
coefficient is the same for the Sun and Ceti the modula- 
tion strength for Ceti can be found by assuming that the 
ram pressure of the wind around Ceti equals that of the 
interstellar medium and therefore is the same as for the Sun: 



pV^ oc 



MV 



(5) 



so 



$ oc r oc V M. (6) 

The mass loss of Ceti has been measured by Gaidos 
(1998) and Gaidos, Giidel & Blake (2000) to M^iceti ~ 
4 • lO-i^Mg/yr. 

Following Shaviv (2003) we can now calculate the en- 
ergy of the GCR influx to an Earth-like planet around 
K^ Ceti relative to the energy of the GCR influx to the Earth: 



(7) 



where Ec is the cutoff energy of GCRs that can actually 
reach the troposphere (~ 12 GeV) and / is the differential 
number flux reaching an Earth-Uke planet around either the 
Sun or K^ Ceti which again is a function of the differential 
number flux of the interstellar medium (Shaviv, 2003): 
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or erticcti/s© ~ 0.1 - i.e. an Earth-like planet around 
Ceti would only receive around 10% of the cosmic ray 
flux that we receive on the Earth and the temperature would 
thus be around 7 degrees warmer than it would have been 
had Ceti not been able to protect its planet more effec- 
tively from OCRs than the Sun. 

3 The effect of Forbush Decreases 

CMEs directed toward the Earth can lead to sudden reduc- 
tions in the influx of OCRs over time scales from hours to 
days known as Forbush decreases. It was shown by Svens- 
mark, Bondo & Svensmark (2009) that large Forbush de- 
creases were followed by reduced levels of aerosols, of 
cloud water content, of liquid water cloud fraction and of 
low IR-detected clouds. 

Svensmark, Bondo & Svensmark (2009) analyzed five 
CMEs found over a time span of 10 years which all re- 
sulted in Forbush decreases associated with an '^lO % de- 
crease in the GCR influx over around a week - this led them 
to note that large Forbush decreases today are too rare to 
have any significant effect on the Earth's climate. On the 
other hand it is not given that large Forbush decreases did 
not have a significant effect on the Earth's climate 4 billion 
years ago when life started to evolve on the Earth, as CMEs 
are expected to have been much more common on the early 
Sun. We therefore analyze if the reduction in the influx of 
OCRs originating from Forbush decreases could be signifi- 
cant compared to the reduction in the influx of OCRs from 
the more effective shielding from the Sun at the time life 
evolved on the Earth. 

Assuming that the CME rate scales liniarly with the flare 
rate (which seems to be the case for the Sun) we can use the 
flare rate of Ceti to provide us with an estimate of the 
CME rate and thus an estimate of how common Forbush 
decreases would be on an Earth-like planet around Ceti. 
The cumulative flare occurrence rate distribution of Ceti 
was calculated by Audard et al. (2000) using 7 days of EUV 
observations from the Extreme Ultraviolet Explorer (Malina 
& Bowyer 1991). In order to compare this cumulative flare 
occurrence rate distribution to the Sun we have analyzed 
the occurrence of solar flares from 1998 to 2007 observed 
in X-ray by the Geostationary Operational Environmental 
Satellite systerrQ (Garcia 1994). 

The flare occurrence rate distribution for the Sun and 

Ceti are shown in Fig. 1 (solid lines) together with power 
law fits to these distributions of the form (Audard et al. 
2000): iV(> E) = kE-"+'^, where fc is a normalization 
factor and a is a constant measuring the hardness of the 
flare distribution. 

In order to compare the two flare occurrence rate distri- 
butions we need to correct them for the fact that the solar 
flares were observed in an energy range roughly 50 times 
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' The Geostationary Operational Environmen- 

tal Satellite system observations were obtained from 
|http://www.ngdc. noaa.gov/stp/SOLAR/ftpsolarflares .html I 
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Fig. 1 Comparison between the cumulative flare occur- 
rence rate distributions for the Sun and Ceti. The flare 
rate for the Sun has been calculated from soft X-ray (1-8 
A) data from The Geostationary Operational Environmen- 
tal Satellite system, and the flare rate for Ceti has been 
calculated from EUV data from the Extreme Ultraviolet Ex- 
plorer (0.01-10 KeV) from Audard et al. (2000). The soHd 
lines show the observations. The dotted lines show power 
law fits to the observations of the Sun and Ceti , respec- 
tivily. The dashed line shows the power law fit to the solar 
observation, but here the flare energies have been multiplied 
by 50 in order to make a reliable comparison to the Ceti 
observations, which have been integrated over a larger en- 
ergy range. It is seen that whereas flares with integrated en- 
ergies around 10^^ erg (the ones that causes Forbush de- 
creases) are rather rare on the Sun, they occure daily on 
Ceti. 



smaller than energy range in which the flares on Ceti 
where observed in; the solar flares were observed in the soft 
X-ray band (1-8 A), while the flares on Ceti were ob- 
served in the EUV band (0.01-10 KeV). We therefore mul- 
tiply the solar flare energies with 50 (Audard et al. 2000). 

All the five Forbush decreases analyzed by Svensmark, 
Bondo & Svensmark (2009) were associated with flares 
with an integrated energy around 10"^^ erg. These Forbush 
decreases generally led to a '^10 % decrease in the GCR 
influx over around a week. In Fig. 1 it is seen that such 
flares happen around once a day on Ceti. This means 
that around 7 Forbush decreases would be present around an 
Earth-like planet around Ceti at any given time. Assum- 
ing that all 7 Forbush decreases will lead to a 10% reduc- 
tion in the GCR influx, an Earth-like planet around Ceti 
would be experiencing an approximate 50% mean reduc- 
tion in the GCR influx from Forbush decreases - i.e around 
half the reduction that is expected to occur due to the more 
effective shielding of GCR around Ceti. 

It is apparent that an Earth-like planet around Ceti 
would be experiencing a reduction in the GCR influx from 
both a more effective shielding from a larger astrosphere 
and from Forbush decreases at the same time. Thus it would 
experience a 90% reduction from more effective shielding 
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from a larger astrosphere and a 50% reduction of the re- 
maining 10% from Forbush decreases. Therefore by adding 
the contribute from Forbush decreases to the contribute 
from more effective shielding from a larger astrosphere 
Eq. 1 predicts an 8 instead of a 7 degree warmer climate. 
This is of course an insignificant difference - i.e. it does not 
change much to remove what is absent. On the other hand 
this study has shown that if the early Sun had not been more 
effective in shielding the Earth from OCRs, the Earth would 
still have experienced a reduced GCR influx due to Forbush 
decreases. 



Unfortunately, we still lack a good asteroseismic esti- 
mate of the age of this star This is unfortunate because 
the best current age estimates of Ceti based on measure- 
ments of the rotation and activity of Ceti come with un- 
certainties of ^ 500 million years, Ceti is therefore an 
obvious target for future ground-based asteroseismic cam- 
paigns - e.g. the first observations by the Stellar Observa- 
tions Network Group (Grundahl et al. 2009). 



4 Conclusion 

Using the young solar twin Ceti as a case study we have 
shown that the reduction in the GCR influx to the Earth 
caused by Forbush decreases had the same order of mag- 
nitude as the reduction caused by more effective shielding 
from a larger heliosphere at the time life evolved on the 
Earth. 

This does not change the conclusion made by Shaviv 
(2003): that the warming associated with a reduced GCR 
influx is enough to significantly compensate for the fainter 
Sun at the time life evolved on the Earth and can explain 
about 1/2 to 2/3's of the temperature increase between now 
and then. Thus the warming associated with a reduced GCR 
influx is enough to solve the faint young Sun paradox. 

An open question is whether the GCRs that are scattered 
away from the CMEs in the Forbush decrease will eventu- 
ally return. There is no evidence from ground-based obser- 
vations of GCRs that the influx increases after a Forbush de- 
crease and it therefore seems secure to assume that the 10% 
of the GCR influx is simply scattered so much away from 
the Earth during a Forbush decrease that it can be consid- 
ered removed from the near-Earth environment. This is most 
likely also the case for an Earth-like planet around Ceti 
as a large part of the CMEs will be magnitudes larger than 
what we have observed on the Sun. 

Ceti is a unique laboratory for understanding the ac- 
tivity of the Sun when life evolved on the Earth. Not only 
for understanding how a reduced GCR influx could affect 
the climate back then, but also for understanding for exam- 
ple how a larger X-ray flux from the larger corona and a 
larger UV flux from a stronger chromosphere could affect 
the evolution of life. 

The larger flare rate of the young Sun has of course 
also had other consequences. Firstly, the largest flares and 
accompanying CMEs would have reduced the amount of 
ozone in the Earth atmosphere, making life on Earth much 
more vulnerable to UV radiation (Schaefer et al. 2000). Sec- 
ondly, large solar flares and accompanying CMEs might 
also have played a more direct role in the evolution of life 
on the Earth by providing an energy source to create organic 
molecules such as the lightning in the Miller-Urey experi- 
ment (Miller & Urey 1959). 
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